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RBC collaborations on QCDOC machines. 
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1. INTRODUCTION AND SIMULATION PARAMETERS 

This work describes results from 2+1 flavour domain wall ensembles with a lattice size of 
16 x 32 recently generated by the RBC and UKQCD collaborations. The analysis was carried out 
using the standard domain wall Dirac operator [p and Pauli-Villars fields with the action introduced 
in The configurations were generated with the Iwasaki gauge action [|J], using the RHMC 
algorithm [Q, ||] with a trajectory length of 1.0, domain wall height of 1.8, j8 = 2.13 and a fifth 
dimension of length 16. Three ensembles were generated with a fixed approximate strange quark 
mass, am s = 0.04, and a light isodoublet with masses am Ui / = 0.01, 0.02 or 0.03. The number of 
trajectories in each ensemble and the number of configurations used in the analysis are shown in 
table HI At present, 24 3 x 64 ensembles with the same parameter values are being generated in 
order to study finite size effects and to compute a wider range of physical quantities. 

In this paper we consider meson correlators with valence quark masses equal to the light 
isodoublet in the sea, the strange quark mass in the sea and the non-degenerate combinations. 
To improve our statistics, correlators were oversampled and averaged into bins of size between 
five and ten, depending on the Monte Carlo time separation between measurements. The binning 
is consistent with the integrated auto-correlation length for the pseudoscalar meson, which was 
calculated to be of order 50 trajectories. Multiple sources per configuration and several different 
types of smearing have also been used to improve the signal. A full correlated analysis was then 
performed with the binned data as input. 



2. RESULTS 



2.1 The residual mass 



The residual mass, which is a measure of the violation of chiral symmetry due to the finite 
fifth dimension was calculated from the ratio of the point-split pseudoscalar density, J$, at the 
middle of the fifth dimension, to the pseudoscalar density, P, built from the fields on the walls [[Tj] 

T. 3 s(J5(y,t)P(x,o)) 

^x,y\P{y,t)P{x,0)) 

A simple unitary linear chiral extrapolation was performed, am ms = A + Bam y, where am f is the 
input valence quark mass and the self-consistent value, am^ es = rrq~gj . was defined to be the residual 
mass in the chiral limit. The self-consistent value obtained in the chiral limit is given in table 
This corresponds to a residual mass of 5 MeV. 
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Table 1: RHMC 2+1 flavour datasets, the value of the residual mass in the chiral limit, am^ es , the bare 
strange quark mass and the unitary pseudoscalar meson mass, amp$. 
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Figure 1: LEFT: Pseudoscalar meson effective masses. Different colours correspond to different quark 
masses, open, closed and square symbols denote different smearings. Horizontal lines show the fit to 



the ground state mass while the curved lines result from inputting the fit parameters into equation (2.2). 



RIGHT: Chiral extrapolation of the pseudoscalar meson mass against average valence quark mass m q = 
(am™ 1 + am q r) /2. The slope shown comes from parameter B in equation (2.4) and corresponds to extrapo 



lation to the chiral limit. Different colours denote different light sea quark masses. It may be seen that the 



sea quark mass dependence, parameter C in equation (2.4), is very mild 



2.2 Light meson masses 

Pseudoscalar and vector meson masses were extracted by performing a simultaneous double 
cosh fit to the form 
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(2.2) 
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to extract both the ground and first excited states. This helped to remove the systematic error in the 
choice of fit range due to the excited state. Figure [I] (left) shows typical effective masses for the 
pseudoscalar mesons. 

We set the scale in three ways. Firstly, the lattice spacing was obtained from m p in the chiral 
limit by performing a simple linear unitary chiral extrapolation of the vector meson mass 



am v =C{am v q al=sea )+D 



(2.3) 



where the quark mass, am q , was defined by am q = a(nif + m^s) with arrif the input quark mass. 
Secondly, ro was determined from the quark-antiquark static potential presented in this confer- 
ence [||], and finally, the method of planes [{7]] was used. All three methods gave a consistent value 
for the lattice spacing of aT x = 1.60(3) GeV, and therefore a box size of ~ 2fm. The lightest 
pseudoscalar meson has a mass of approximately 390 MeV. 

The non-degenerate pseudoscalar meson masses were chirally extrapolated by fitting both the 
valence and sea quark mass dependence using the form 



{amp S ) 2 =A + B (am q f + am q f) + Cam 



sea 

q ■ 



(2.4) 
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Figure fright) shows the pseudoscalar meson mass squared plotted against the average valence 
quark mass, am^ = (am™ 1 + am™ l )/2. The value of A was found to be consistent with zero in good 
agreement with the residual mass behaving as an additive renormalisation to the input quark mass, 
nif. The normal quark mass, am n , was obtained from the physical pion mass by setting am n = 
am™ 1 = am™ 1 = am^ a in equation (^4j). Similarly, the strange quark mass, am s , was obtained from 
the physical kaon mass by substituting am n = am™ 1 = am s q ea and am s = am™ 1 . In both cases the 
scale was taken to be 1.60(3) GeV. The value obtained for the bare strange quark mass (see table p]) 
is in agreement with the input value of 0.04 + am?, 



'res- 



2.3 Non-perturbative renormalisation for Bk 

Bk measures the QCD correction to the weak mixing between K° = (sd) and ^° = (sd) which 
is an important constraint on the CKM matrix. It is evaluated from the dimensionless ratio of the 
matrix element of the AS = 2 operator, 

@b k = 0vv+aa = (sr^d)(sY^d) + (sY 5 Y^d)(sY 5 Y^d) , (2.5) 

to its vacuum saturation approximation 

<jft|<y,.|g°) 

Chiral symmetry breaking may cause mixing of the operator with four other four-quark oper- 
ators 

0VV-AA = [sYiid){sYiid) - [sY5Yiid){sY5YiJid), (2.7) 

& SS±PP = (sd) (sd) ± (sYsd) (sjsd) , (2.8) 

Gtt = (sa^ v d)(sa^ v d) , (2.9) 

however, the exponentially accurate chiral symmetry afforded by the domain wall fermion for- 



mulation should strongly suppress this mixing. In the presence of chiral symmetry breaking the 
renormalised B-parameter, B K en ^ , is in principle an admixture of the five lattice fi-parameters 

R (ren) _ y ,„% R (latt) , y „ /y 2 R (latt) / 91m 

ti K — LvV+AA,VV+AAI ^A' Br + l_, L VV+AA,il ^A ' "i > U-AU) 



for i € {VV —AA,SS ±PP,TT}. It may be observed from figure^ (left) that the magnitude of the 
B-parameters for the wrong chirality operators are large in comparison with By V+AA , the matrix 
element of interest, and hence the sizes of the mixing coefficients, Zw+aaj, are crucial. 

The renormalisation coefficients are evaluated using the Rome-Southampton method []|] of 
non-perturbative renormalisation. In Landau gauge, the amputated n-point correlation function of 

( n ) ltltt 

the operator of interest, TV , is constructed by applying, at fixed external quark momenta and 
zero quark mass, the condition 

r ^ latt (^;a- 1 )]=P,[rS; )tree ], (2.11) 
where represents the application of a particular spin-color projection and a subsequent trace. 



Z q n/2 ZijP k 
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Figure 2: LEFT: B-parameters for 



'vv 



-AA, &VV-AA, &SS+PP, &SS-PP an d 



Piatt 

Gjt as a function of input 



RIGHT: Renormalisation factors Z„ 2 Zy 



as a function of p^. This demonstrates 



quark mass, mf. Kium; K.eiiormaiisauoii laciors z,^ ~^vV+AA.j 
that the mixing matrix is, in practice, block diagonal with negligible chirality mixing. For each momentum 
the data has been linearly extrapolated to the chiral limit (m q — » 0). 



The set of momenta used to calculate the ratio of renormalisation factors Zq(aj=2) /Z\ is defined 

by 

(2k 2tl 2% 2k \ 

Piatt = —n x , —n y , —n z , —n t (2.12) 
\ L, x Ly L z L t J 

where L x = L y = L z = 16 and L t = 32, using 1359 combinations of (n x ,n y ,n z ,n t ) with —2 < 
n x ,n y ,n z < 2 and — 5 < n t < 5. These are then averaged into 29 combinations of equal pf M . The 
number of gauge configurations used for the non-perturbative renormalisation is given in table [T| 

Results for the elements of the top row of the mixing matrix, Aw+aaj, ^q^vv+AA p versus P\ a tt 
are shown in figure || (right). It may be seen from the plot that, in the window of momenta for which 
contributions from both hadronic effects (low momenta) and contributions from discretisation ef- 
fects (high momenta) are small, the resulting mixing coefficients are small enough that the matrix 
becomes block diagonal, simplifying the inverse, and hence the wrong chirality operators can be 
safely neglected. Therefore, we calculate the Bk renormalisation factor as z^/ MOM = Zw+aa/Z,^. 
In order to obtain a value for Bk in the MS scheme at some scale a continuum running/matching 
calculation was performed following the techniques in [g]. A preliminary value for Zg K {MS, 2GeV) 
of 0.90(2) was obtained. 



3. SUMMARY 

Recent 2+1 flavour simulations using domain wall fermions, the Iwasaki gauge action, three 
different light isodoublet quark masses and a volume of 16 3 x 32 with a fifth dimension size of 
16 have been performed by the RBC and UKQCD collaborations. The chiral symmetry breaking, 
parameterised by the residual mass, is consistent with a small additive quark mass renormalisation 
of 5 MeV, which is at an acceptable level with a fifth dimension size of 16. Reasonable signals 
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have been obtained for the pseudoscalar and vector meson masses. These ensembles have a lattice 
spacing of 1.60(3) GeV and hence the lightest pseudoscalar meson has a mass of approximately 390 
MeV. We find that the bare strange quark mass evaluated from the pseudoscalar chiral extrapolation 
is in good agreement with the input strange quark mass. 

The improved chiral symmetry afforded by the domain wall fermion formulation leads to the 
suppression of mixing with the wrong chirality operators in the calculation of Bk- This greatly sim- 
plifies this calculation and should allow for an accurate determination of the matrix element [[R]] . 
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